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ABSTRACT 

The  problem  of  determining  the  optimum  allocation  of  aircraft  to 
an  airstrike  against  a  transportation  network  is  investigated.  The 
damage  function  is  assumed  to  be  exponential.  A  solution procedure^  is 
developed  utilizing  dynamic  programming  and:  integer  solutions,  aree found, 
The  number  of  aircraft  to  be  assigned  to  the:  airstrike  is:  considered  a 
decision  variable.  A  sensitivity  analysis  is  run  to  determine- the 
optimum  value  for  this  variable. 
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I.      INTRODUCTION 

A.     OBJECTIVE 

The  purpose  of  this  paper  is  to  present  a  procedure  for  determining 
tiree  optimal   allocation  of  aircraft  to  a  single  airstrike  against  a 
transportation  network.     This  allocation  problem  is  solved  by  dynamic 
programming  and: a  Fortran-coded  version  of  the  program  is  included  in 
ttre  paper . 

ft.     GENERAL 

Sustained  ground  operations  require  a  military  force  to  have  some 
HEffiBTs:  of"  resupply.  This  resupply  capability  is  partially  dependent 
upon  a  land  transportation  system.  The  level  of  resupply  effort  required 
depends  upon  what  type  of  forces  are  being  supported.  Guerrilla  forces 
enjoying  local  support  require  less  resupply  capability  in  terms  of 
pounds  per  man  per  day  than  would  a  conventional  army,  but  a  greater 
percentage  of  this  capability  depends  upon  land  transportation  networks. 

Any  reduction  in  the  resupply  capability  of  a  military  force  will 
reduce  its  combat  effectiveness.  Tactical  air  interdiction  has  been 
used  extensively  by  the  Armed  Forces  of  the  United  States  against  its 
opponents  in  Southeast  Asia  to  accomplish  this  reduction. 

There  are  at  least  three  alternative  means  of  using  tactical  air  to 
reduce  the  resupply  capability  of  an  enemy.  Aircraft  may  be  assigned 
to  attack  sources  of  supply  to  destroy  war  material  before  it  enters 
the  transportation  system  and/or  to  disrupt  its  production;  aircraft 
caai  destroy  war  material  as  it  moves  in  the  transportation  system;  and 
finally  aircraft  can  attempt  to  reduce  the  resupply  capacity  of  the 
transportation  system  itself  by  destroying  bridges,  roads,  railroads, 


et  cetera.  Conventional  wisdom  argues  that  the  first  course?  o~f  action 
is  the  most  effective  form  of  interdiction.  Unfortunately  for  military 
planners,  political  considerations  may  rule  out  this  alternative.  This 
paper  will  focus  on  the  last  of  these  options,  the  reduction  in  capacity 
of  the  transportation  system  itself. 

C.  BACKGROUND 

Considerable  effort  has  been  devoted  to  the  interdiction  problem. 
In  particular  two  recent  papers  provided  the  background  for  the  approach 
to  the  problem  developed  in  this  paper.  McMasters  and  Mustin  [1] 
developed  an  algorithm  that  determines  which  arcs  of  a  transportation 
network  should  be  attacked  and  at  what  level  of  effort  giverr  a:  limited 
availability  of  resources.  In  this  formulation  of  the  problem  the 
relationship  between  arc  capacity  and  resource  allocation  (damage 
function)  was  assumed  to  be  linear.  The  algorithm  presented  is  based 
upon  the  max- flow  min-cut  theorem  of  Ford  and  Fulkerson  [2]  and  the 
relationship  between  a  primal  network  and  its  topological  dual. 

Nugent  [3]  investigated  the  same  problem  under  the  assumption  of  an 
exponential  damage  function  which  exhibits  diminishing  marginal  returns. 
An  algorithm  was  developed  that  finds  a  non-integer  solution  to  the 
problem. 

In  this  paper  the  transportation  system  will  have  the  same  network 
formulation  as  in  Refs.  1  and  3.  The  problem  will  be  formulated  differ- 
ently and  dynamic  programming  will  be  used  to  provide  integer  solutions. 

D.  INTERDICTION  PROBLEM 

It  will  be  assumed  that,  given  unlimited  aircraft  availability,  the 
assignment  of  aircraft  to  an  airstrike  would  reach  a  point  beyond  which 
it  would  become  uneconomical  to  assign  further  aircraft.  In  a  problem 


with  constraints  on  aircraft  availability  this  point  might  or  might  not 
occur  before  all  available  aircraft  were  assigned.  For  this  reason, 
the  objective  of  an  operations  officer  planning  an  atrstrike  against  a 
transportation  network  is  not  merely  to  minimize  network  capacity 
subject  to  aircraft  availability,  but  to  minimize  the  capacity  subject 
to  aircraft  availability  and  the  additional  consideration  that  the  cost 
of  any  incremental  assignment  of  aircraft  to  the  strike  ts  exceeded  by 
the  benefit  resulting  from  that  assignment. 

To  accomplish  the  objective  the  strike  planner  must  have  information 
on  the  availability  and  cost  of  assignment  of  aircraft,  detailed  infor- 
mation must  be  available  concerning  the  transportation  network  including 
the  upper  and  lower  bounds  on  the  capacity  of  each  arc  and  its  vulnera- 
bility to  attack.  The  planner  must  also  know  the  benefit  to  attribute 
to  a  reduction  in  resupply  capability.  With  this  information  and  using 
the  procedure  that  will  be  outlined  the  planner  can  determine:  how 
many  aircraft  to  assign^  to  the  airstrike;  which  arcs  in  the  network 
should  be  attacked;  how  many  aircraft  to  assign  to  arcs  that  will  be 
attacked;  and  the  capacity  of  the  network  after  the  airstrike. 


II.  THE  MODEL 

A.  NETWORK  DESCRIPTION 

The  transportation  system  under  consideration  is  represented  by  a 
planar  connected  graph  of  nodes  and  undirected  capacitated  arcs.  Arcs 
represent  road  segments  and  nodes  represent  either  a  read  irrtersection 
or  any  other  point  where  it  is  necessary  to  distinguish  between  road 
characteristics  on  either  side  of  the  node.  Three  constants  are- asso- 
ciated with  each  arc  representing  the  upper  and  lower  bounds  on  arc 
capacity  and  the  arc's  vulnerability  parameter. 

It  is  assumed  that  the  network  has  one  source  node  from  which  flow 
originates  and  one  sink  node  at  which  flow  terminates.  If  the  trans- 
portation system  being  modeled  has  more  than  one  originating  point  or 
terminating  point  this  may  be  handled  by  creating  a  super-source  and/or 
sink  with  artificial  arcs  connecting  these  super-nodes  to  sources  and 
sinks  as  needed.  These  artificial  arcs  may  not  be  attacked  and  their 
capacities  are  unbounded.  The  arc  between  nodes  i  and  j  is  represented 
by  (i,j).  Nodes  are  numbered  from  1  to  n  with  1  corresponding  to  the 
source  and  n  the  sink.  With  the  exception  of  the  source  and  the  sink, 
flow  conservation  is  assumed  to  hold.  That  is,  flow  out  of  node  i 
equals  flow  into  node  i. 

The  flow  in  arc  (i,j)  is  designated  as  x. .  if  it  is  from  node  i  to  j 
and  x-.j  if  it  is  from  node  j  to  i .  This  avoids  the  necessity  of  defining 
negative  flows.  Flow  is  assumed  to  be  from  the  source  to  the  sink 
although  it  may  be  in  either  direction  in  the  intermediate  arcs.  The 
model  as  formulated  considers  only  flows  of  a  single  commodity,  tons 


of  resupply  per  day,  and  the  value  of  one  unit  of  flow  is  assumed  to 
be  the  same  for  all  arcs. 

Capacities  on  arcs  represent  bounds  on  flow  in  either  direction. 
The  capacity  on  arc  (i,j)  is  given  by  m. .  and  is  assumed  to  be  the  same 
in  both  directions.  The  flow  in  arc  (i,j)  is  restricted  by 

0  1  xlj  <  mi  j  ' 

The  upper  and  lower  bounds  on  the  capacity  of  arc  (i,j)  are  repre- 
sented by  u- •  and  1 . .  where 

0  <  1  ,•  n-  <  m.-  n-   <  u .  .    . 
—    1J  —    U  —    ij 


th 


The  vulnerable  portion  of  an  arc's  capacity  is  designated  w.jj  wi 

Wi  i     =     u . .  -  1 . .   . 

The  amount  of  resource  allocated  to  interdict  arc  (i,j)  is  denoted  by 
Lsj.  The  relationship  between  the  capacity  of  arc  (i,j)  and  the  level 
of  resource  assigned  to  its  interdiction  is  defined  as  the  damage 
function  of  arc  (i,j)  and  is  given  by 

mtj(ki;j)  =  lia.  +Wij  expf-byky)  . 

In  the  above  damage  function  the  parameter  b-.  is  a  measure  of  the 

i  j 

vulnerability  of  arc  (i,j).  Larger  values  of  b--  result  in  greater 

reductions  in  capacity  for  fixed  values  of  I,--.-,  w--  and  k . .  and  hence 

•j   ij     ij 

imply  greater  vulnerability.  If  bjj  =  0  then  m^-Ck.  .)  =  un-.-  for 
all  possible  values  of  k^j  and  no  capacity  reduction  is  possible.  With 
this  damage  function  if  no  aircraft  are  assigned  to  (i,j),  its  capacity 
will  be  uy  and  in  the  limit  as  the  number  of  aircraft  assigned  to  (i,j) 


becomes  infinite  the  capacity  approaches  1^.     This  lower  bound  will 
be  referred  to  as  arc  capacity  after  unlimited  interdiction. 


B.  DETERMINATION  OF  NETWORK  CAPACITY 

The  opposition  is  assumed  to  have  the  means  to  determine  how  to 
maximize  the  flow  in  the  transportation  network..     Let  the:  capacity  of 
the  network  be  defined  as  this  maximal   flow.     The  determination  of 
maximum  flow  is  the  well-known  maximal   flow  problem  and  may  be  found 
using  the  max-flow  labeling  algorithm  based  upon  the  max-flow  min-cut 
theorem  of  Ford  and  Fulkerson  [2].     Ford  and  Fulkerson's  theorem  states 
that  the  maximum  flow  possible  in  a  network  is  equal   to  the  value  of 
the  minimal   cut  set.     In  this  paper  the  value  of  a  cut  set  will  be 
referred  to  as  its  capacity. 

C.  ENUMERATION  OF  CUT  SETS 

The  network  capacity  has  been  defined  to  be  equal  to  the  maximum 
flow  possible  in  the  network.  As  discussed,  this  maximum  flow  is  equal 
to  the  value  of  the  minimum  cut  set.  Therefore,  the  problem  of  mini- 
mizing this  capacity  is  equivalent  to  minimizing  the  capacity  of  some 
cut  set.  It  is  obvious  that  aircraft  will  be  allocated  to  only  one  cut 
set  since  if  this  were  not  the  case  all  aircraft  could  have  been  assigned 
to  the  cut  set  that  was  minimal  after  the  first  allocation  with  a 
resulting  decrease  in  network  capacity. 

The  complicating  factor  is  that  there  is  no  easy  way  to  find  out 
which  cut  set  should  be  selected  for  attack.  To  solve  the  problem  it 
is  necessary  to  have  some  means  of  identifying  cut  sets.   In  addition, 
it  is  desirable  to  be  able  to  identify  these  cut  sets  in  order  of 
increasing  capacity  after  unlimited  interdiction  since  once  a  cut  set 
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is  found  whose  capacity  after  unlimited  interdiction  is  greater  than  or 
equal  to  network  capacity  before  interdiction  no  more  cut  sets  need  be 
identified.  The  network  capacity  before  interdiction  represents  an 
upper  bound  on  network  capacity.   Define  S^  as  the  cut  set  with  the  itn 
smallest  capacity  after  unlimited  interdiction.  The  set  of  S.  whose 
capacities  after  unlimited  interdiction  is  less  than  the  upper  bound  on 
network  capacity  will  be  denoted  by  S. 

The  method  by  which  cut  sets  are  identified  makes  use  of  the  topo- 
logical dual  of  a  network.  Arcs  have  lengths  rather  than  capacities  in 
the  dual  network.  The  cut  sets  in  the  primal  network  have  a  one-to-one 
correspondence  with  the  loop! ess  paths  in  the  dual.  The  problem  of 
finding  the  shortest  path  from  the  dual  source  to  the  dual  sink  corre- 
sponds to  the  primal  problem  of  finding  the  minimum  cut  set.  The  length 
of  the  dual  shortest  path  equals  the  primal  capacity. 

The  topological  dual  of  a  given  primal  network  is  constructed  as 
f ol 1 ows : 

(1)  Connect  the  source  and  the  sink  of  the  primal  with  an  artificial 
arc.  Call  the  result  the  modified  primal. 

(2)  Place  a  node  in  the  area  surrounding  the  modified  primal  (external 
face)  and  one  in  each  face  formed  by  the  arcs  of  the  modified  primal. 
Let  the  dual  source  be  the  node  in  the  external  face  and  the  dual  sink 
be  the  node  in  the  face  involving  the  artificial  arc. 

(3)  For  each  arc  in  the  primal  (except  the  artificial  arc)  construct 
a  dual  arc  that  intersects  it  and  joins  the  two  nodes  in  the  faces 
adjacent  to  it. 

(4)  Assign  each  dual  arc  a  length  equal  to  the  capacity  of  the  primal 
arc  it  intersects. 
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Once  the  dual  network  has  been  developed,,  the  shortest  path  through 
the  dual  before  interdiction  is  found.  This  path  is  determined  using 
the  upper  bounds  on  primal  capacities  as  lengths  of  arcs  in  the  dual. 
The  length  of  this  path  represents  network  capacity  before  interdiction. 
Any  shortest  path  algorithm  may  be  used  for  this  determination.  Dreyfus 
[4J  evaluated  several  of  these  algorithms  concluding  that  the  procedure 
developed  by  Dijkstra  is  the  mast  efficient.  Next  the  lengths:  of-  the: 
dual  arcs  are  changed  to  correspond  to  the  lower  bounds  on  primal  arc 
capacities.  The  lengths  of  the  dual  paths  now  represent  the  capacities 
of  the  corresponding  primal  cut  sets  after  unlimited  interdiction. 
Paths  with  loops  need  not  be  considered  since  they  correspond:  to  primal 
cut  sets  that  either  include  more  arcs  than  necessary  to  sever  the 
network  or  contain  some  arc  more  than  once.  The  dual  paths  are  identi- 
fied in  order  of  increasing  length  by  means  of  an  ntn  shortest  path 
algorithm.  Clarke,  Krikorian  and  Rausen  [5]  developed  an  algorithm  for 
determining  the  n  best  loopless  paths,  but  it  is  difficult  to  apply. 
Pollack  [6]  in  an  unpublished  paper  presented  an  algorithm  which  succes- 
sively develops  the  best  loopless  paths  using  extensions  of  shortest 
path  algorithms.  This  procedure  is  less  complex  than  that  of  Clarke, 
Krikorian  and  Rausen  and  appears  to  be  more  efficient.  It  should  be 
noted  that  depending  on  the  number  of  elements  in  S  and  the  total 
number  of  paths  in  the  dual,  the  most  efficient  means  of  developing  S 
may  be  to  enumerate  all  paths  through  the  dual  and  then  compare  lengths. 


12 


III.  ANALYSIS  OF  THE  MODEL 

A.  MATHEMATICAL  FORMULATION 

The  problem,  as  outlined  previously,  is  to  find  that  allocation  of 
aircraft  to  an  airstrike  against  a  trans partat ion  network  which  will 
minimize  the  capacity  of  that  network.  This  mtntmizatxon  is:  accomplished 
subject  to  a  constraint  on  aircraft  availability  and  the.  consideration 
that  the  incremental  benefit  of  assigning  aircraft  must  exceed  the 
Incremental  cost  of  that  assignment.  If  net  benefit  is  defined  to  be 
the  difference  between  the  total  benefit  derived  from  ther  airstrike  and 
the  total  cost  of  aircraft  assignment  the  problem  may  b&  restated  as 
follows:  maximize  net  benefit  subject  to  aircraft  availability. 

Let  K  represent  the  total  number  of  aircraft  available  for  assign- 
ment to  the  airstrike  and  let  K*  be  the  number  of  aircraft  that  have 
been  assigned  to  the  airstrike.  Then  for  any  choice  of  K*  the  problem 
may  be  stated  mathematically  as 

min   [cut  set  capacity  after  optimal  interdiction] 


Vs 


or 


min       [  min        z  (1..  +w..  exp{-b,.k..} 

SjeS  (1,j)eSi     1J         ]J  ^   *J 


)] 


subject  to  E  k, .  <  K* 

CtJks1    1J~ 

k.jj       positive  integer 
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The  structure  of  this  problem  will  allow  the  devel  apment.  of  an 
efficient  solution  procedure.  Note  that  with  respect  to  a  particular 
cut  set  the  objective  is  to  minimize  its  capacity.  Since  the  cut  set 
capacity  is  the  sum  of  functions  that  are  convex  in  k-j-,  this  capacity 
is  a  convex  function  and  is  therefore  unimodal  with  respect  to  minimiza- 
tion. The  overall  objective  function  is  the  minimum  of  a  set.  of- convex 
functions  and  is  neither  concave  nor  convex.  This:  together  with- the 
problem  of  not  knowing  which  cut  set  is  going  to  be  attacked:  requires 
that  each  cut  set  in  S  be  the  subject  of  a  minimization  problem. 

For  a  particular  cut  set,  S-j ,  the  problem  is 

min   E     (1*1  +  *..-  exp[-bf,k.  ]) 
(i,j)eS.   1J    1J      TJ  TJ 

subject  to    E      k.-j  <  K 

k-.  positive  integer  . 

The  term   E  1.,.  is  constant  and  may  be  deleted  during  the  minimization 

(i,j)  1J 
and  then  added  back  to  give  the  solution  in  terms  of  capacity.  This 

problem  will  be  solved  by  means  of  dynamic  programming.  Each  arc  in 

the  cut  set  under  consideration  will  be  represented  by  a  stage  in  the 

dynamic  program. 

Let  the  number  of  arcs  be  n  and  resubscript  each  arc  (i,j)  and  its 

associated  parameters  in  any  order  with  the  single  subscript  i  running 

from  one  to  n.  The  decision  variable  for  stage  i  is  k-  and  the  return 

function  for  stage  i  is  given  by 


^  =  wi  expC-b^)  . 
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The  state  variable  for  stage  i  will  be  denoted  try  x^.  an±  represents:  the 
remaining  resource  availability  at  stage  i.  The  problem  may  be  restated 
as  finding  fnUn)  where 

n 

fn(Xn}     =    n1'?  «Zi   r1(xtJ 

Olknlxn     1=1 

subject  to        x.  -r  =  x.  -  k.     . 

F  (xR)  is  the  optimal  return  from  stages  n,n-l  ,.  ..,1  given  xn  units  of 
resource.  The  above  problem  may  be  solved  by  dynamic  programming  since 
fn(xn)  can  be  decomposed  into  a  series  of  single  variable;  problems. 
Nemhauser  [7]  shows  that  problems  wfth  addtttve  s±age:  returns:  may  always 
be  decomposed.  Therefore,  the  following  recursive  relationship  is  valid 

fn(xn)  =  min     [rn(kn)  +  f^U^)] 
0<^kp<xn  ,  n>l 

and 

f](x-,)  =  min     ri(xl)  • 

0<k]<X] 

The  value  of  xn_-j  is  given  by 

xn_l  =  ^n^n'^n'  =  xn  "  ^n 

where  tR  is  the  transformation  which  gives  the  relationship  between  the 
amount  of  resource  remaining  after  stage  n  given  that  xn  was  available 
before  stage  n  and  kn  was  utilized  at  stage  n. 

The  dynamic  program  is  solved  by  starting  at  stage  one  and  working 
to  stage  n  solving  a  series  of  single  variable  minimizations.  These 
minimizations  are  facilitated  by  the  convexity  of  ttre  individual  stage 
returns.  Nemhauser  [7]  provides  a  proof  of  the  fact  that  in  the 
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minimization  of  additive  stage  returns  the  convexity  of  each  stage 
return  ensures  that  fA*j)   is  a  convex  function  of  xi .  This  means  that 
each  single  variable  optimization  performed  in  the  dynamic  program  is 
of  a  unimodal  function  and  permits  the  use  of  Fibonacci  search  to  find 
the  optimal  values  of  the  decision  variables.  An  application  of- this 
technique  is  found  in  Ref.  7. 

After  the  optimal  allocation  of  aircraft  within-  each  cutset:  in  S 
is  found  for  a  given  K*,  their  capacities  are  compared.  The  cutset 
with  the  minimum  capacity  is  the  one  that  would  be  attacked  if- K*"  air- 
craft were  to  be  assigned  to  the  airstrike.  The  capacity  of  this' minimal 
cut  set  is  by  definition  the  network  capacity  and  this:  capacity  will  be  a 
strictly  decreasing  function  of  K*.  The  remaining  problem  is  to  deter- 
mine how  many  aircraft  to  assign  to  the  strike  in  order  to  maximize 
net  benefit.  To  make  this  determination  it  is  necessary  to  know  the 
cost  of  allocating  aircraft  to  the  strike.  This  will  be  assumed  to  be 
a  constant  C  dollars  per  aircraft.  The  benefit  derived  from  network 
capacity  reduction  must  also  be  known.  It  will  be  assumed  to  be  a 
constant  D  dollars  per  unit  capacity  reduction. 

With  the  above  information  the  problem  of  determining  how  many 
aircraft  to  allocate  may  be  determined  by  comparing  the  incremental 
cost  of  assigning  aircraft  to  the  benefit  resulting  from  that  assign- 
ment. To  make  this  comparison  it  is  necessary  to  define  the  benefit 
resulting  from  the  assignment  of  a  single  aircraft.  This  will  be  de- 
fined as  the  product  of  the  benefit  per  unit  capacity  reduction  (D) 
and  the  amount  of  capacity  reduction  that  can  be  achieved  by  that 
aircraft. 
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The  amount  of  capacity  reduction  that  can  be  achieved  by  one 
additional  aircraft  is  a  function  of  the  number  already  assigned  and 
will  be  denoted  as  6(K*).  A  simple  decision  rule  is  to  assign  aircraft 
K*  =  1,2,...  until  a  point  is  reached  where  benefit  from  the  last 
aircraft  assigned  does  not  exceed  the  cost  of  assignment.  At  this,  point 

D  *  5(K*)  <_    C 
or 

6(K*)  <  C/D 

and  the  optimal  allocation  of  aircraft  is  K*-l  .  If  6(K*)  >  C/Q 
for  all  K*  the  optimal  allocation  is  K  under  this  rule.  There  would  be 
no  problems  with  this  decision  rule  if  <5(K*)  were  a  non-i rrcreasi ng- 
function  of  K*.  In  this  case  once  a  K*  was  found  such  that  6(K*)  <_  C/D 
the  cost  of  any  further  assignment  of  aircraft  would  exceed  its  benefit. 

If  network  capacity  after  optimal  interdiction  were  determined  by 
only  one  cut  set  for  all  values  of  K*  then  6(K*)  would  be  non-increasing. 
This  is  not  the  case.  In  general  as  K*  ranges  from  0  to  K  different  cut 
sets  are  minimal  (see  Figure  1).  At  K*  =  0  the  cut  set  that  determines 
network  capacity  is  by  definition  the  one  that  is  minimal  before  any 
interdiction  takes  place.  Unless  this  cut  set  is  also  minimal  after 
unlimited  interdiction,  at  some  point  another  cut  set  must  determine 
network  capacity.  This  crossover  may,  of  course,  occur  after  assignment 
of  all  available  aircraft.  These  points  where  a  change  in  the  constrain- 
ing cut  set  occurs  represent  points  where  6(K*)  increases  with  respect  to 
K*.  Therefore,  there  is  no  guarantee  that  stopping  when  6(K*)  <_  C/D 
for  the  first  time  is  optimal.   If  at  some  point  after  further  assign- 
ment  of  aircraft  is  made  <$(K*)  again  exceeds  C/D,  it  may  be  that  further 
assignment  of  aircraft  would  have  resulted  in  benefits  outweighing  costs. 
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Figure  1.  Network  Capacity 

The  problem  of  determining  the  optimal  K*  will  be  handled  as  follows: 

(1)  Find  the  first  value  of  K*  for  which  5(K*)  <_  C/D.  Subtract  one 
aircraft  and  let  the  resulting  value  of  K*  be  K-|*.   If  K*  =  K  before 
ICj*  is  found  then  the  optimal  allocation  of  aircraft  is  K. 

(2)  Check  to  see  if  6(K*)  >  C/D  for  any  values  of  K*  >  K}*.     If  not 
go  to  step  (4).  If  so  find  the  next  value  of  K*  for  which  6(K*)  <_  C/D. 
Let  this  number  minus  one  be  l<2*. 

(3)  Continue  in  this  manner  to  identify  the  K*  at  which  6(K*)  becomes 
<  C/D  after  there  has  been  an  intervening  value  of  K*  such  that 
6(K*)  >  C/D.  Subtracting  one  aircraft  each  time,  label  the  resulting 
values  K3*,K4*,  ...,  K^*.  If  6(K)  >  C/D  let  Kn*  =  K.  Let  Kq*  be 
defined  as  0. 

(4)  Starting  with  i  =  1  and  continuing  until  i  =  n,  check  whether  or 
not  the  cost  to  reach  K-j  from  K^-j  is  exceeded  by  the  benefit.  If  it  is 
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let  K*   =  K,-*  ,  increment  i  by  one,  and  go  to  the  beginning  of  step 
opt    ' 

(4).  If  it  is  not,  go  to  step  (5). 

(5)  Starting  with  1  =  1  and  continuing  until  1  =  n-i  check  whether  the 
cost  to  reach  K^+-j  from  K^_ -j  is  exceeded  by  the  benefit.  If  it  is,  let 
K*  .  =  K*i+1  ,  let  i  =  i+1+1  and  go  to  step  (4).  If  not,  increment  1 
by  one  and  go  to  the  beginning  of  step  (5). 

At  the  end  of  this  procedure  K*0  t  will  be  the  optimal  numtrer  of  air- 
craft to  assign  to  the  airstrike  and  the  problem  will  be  solved. 


B.  STEPWISE  SOLUTION  PROCEDURE 

(1)  Formulate  the  topological  dual  of  the  transportation  network.  Find 
the  shortest  path  through  the  dual  before  interdiction.  "Hits  represents 
an  upper  bound  on  network  capacity. 

(2)  Use  Pollack's  algorithm  [6]  to  identify  the  first,  second,  third, 

shortest  paths  through  the  dual  using  the  lower  bounds.  Continue 

identifying  paths  until  one  is  found  whose  length  exceeds  the  previously 
found  upper  bound  on  network  capacity.  Let  the  primal  cut  sets  corre- 
sponding to  these  paths  be  denoted  as  set  S. 

(3)  For  each  cut  set  that  is  an  element  of  S,  use  dynamic  programming 
to  find  the  optimal  allocation  of  aircraft  and  the  resulting  capacity 
for  K*  equal  to  1,2,...,K. 

(4)  For  each  value  of  K*  find  the  network  capacity  by  taking  the  minimum 
of  the  capacities  of  the  elements  of  S. 

(5)  Construct  the  function  6(K*)  and  determine  Kj*,  K2*,...,Kn*. 

(6)  Using  the  procedure  previously  outlined  determine  which  K^*  is 
optimal . 
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C.  SAMPLE  PROBLEM 

The  diagram  in  Figure  2  represents  a  hypothetical  transportation 
network.  The  three  numbers  associated  with  each  arc  are  b-jj,  1-jj, 
and  Uij. 


(.20,475,5 


(.24,0,130) 


(.18,320, 


00,400) 


Figure  2.  An  Example  Network 
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Figure  3  shows  how  the  topological  dual  of  the  transportation 
network  is  constructed. 


„—-*•-- 


—  ■3*-- 


Figure  3.  Construction  of  the  Dual 
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Figure  4  shows  the  topological  dual  after  the  data  for  each  arc 
has  been  transferred  from  the  primal.  In  the  dual  u ■ •  and  1^-  represent 
bounds  on  arc  length. 


(.20,475,550)  /   (.22 


(.15,280,405; 


200,280)     (.30, 270, 3^0)   (.>£,320,440 


(.24,0,1350    (.38 


160,280)      (.25,200 
(.04,400,700) 


290) 
(J&8, 190,400) 


Figure  4.  The  Topological  Dual 


To  simplify  notation,  paths  through  the  dual  will  be  designated  by  the 
nodes  over  which  they  pass.  The  shortest  path  through  the  dual  before 
interdiction  is  1,2,5,9  with  a  length  of  1395.  This  gives  an  upper  bound 
on  network  capacity.  Table  I  lists  all  loopless  paths  through  the  dual 
in  order  of  length  after  unlimited  interdiction.  It  should  be  noted 
that  the  length  of  path  number  11,  the  11th  shortest  path  after  unlimited 
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TABLE   I.      DUAL  PATHS 

PATH  NODES 

1  1,4,7,6,9  780 

2  1,4,7,8,9  890 

3  1,4,3,2,5,9  9350: 

4  1,2,5,9  T075 

5  1,4,3,2,5,6,9  1090 

6  1,3,2,5,9  1T50 

7  1,4,7,6,5,9  1190 

8  1,2,5,6,9  mm 

9  1,3,2,5,6,9  1280 

10  1,3,4,7,6,9  1290 

11  1,3,4,7,8,9  1400 

12  1,4,3,2,5,6,7,8,9  1600 

13  1,2,3,4,7,6,9  1615 

14  1,3,4,7,6,5,9  1700 

15  1,2,5,6,7,8,9  1715 

16  1,2,3,4,7,8,9  1725 

17  1,3,2,5,6,7,8,9  1790 

18  1,2,3,4,7,6,5,9  2025 
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interdiction,  exceeds  the  upper  bound  on  network  capacity.  Therefore 
the  cut  sets  comprising  set  S  correspond  to  paths  1  through  10. 

It  is  assumed  for  purposes  of  this  example  that  there  are  100  air- 
craft available  for  assignment  at  a  cost  of  30,000  dollars  for  each 
aircraft  assigned.  It  is  further  assumed  that  the  benefit  derived  from 
a  reduction  of  one  ton  per  day  in  network  capacity  is  7..5QQ  dollars. 

The  dynamic  program  for  each  S^  that  is  constraining  including  a 
sensitivity  analysis  on  K*  is  contained  in  the  computer  output-  A 
graph  of  the  resulting  network  capacity  is  given  by  Figure  5.  For  K* 
in  the  range  1  through  25  cut  set  4  determines  network  capacity,  for 
K*  in  the  range  26  through  57  cut  set  3  is  constraining,  and  for  K*  from 
58  to  100  cut  set  1  is  minimal. 

From  the  given  values  of  C  and  D.  30,000  and  7,500  respectively, 
the  points  of  interest  are  those  at  which  6(K*)  becomes  <_  C/D  =  4.  This 
occurs  for  the  first  time  when  K*  =  42.  Therefore,  Kj*  =  41.  At  K*'  =  58 
6(K*)  again  exceeds  4  so  it  is  necessary  to  search  for  another  point 
where  s(K*)  <  4.  This  next  occurs  at  K*  =  62  and  K2*  is  61.  Since 
6(K*)  does  not  exceed  4  for  any  K*  >  62,  Kn*  =  K2*. 

It  is  obvious  that  the  benefit  to  get  to  K-j*  exceeded  the  cost 
since  K-j*  was  the  first  point  at  which  the  allocation  of  another  aircraft 
did  not  produce  benefits  exceeding  costs.  However,  it  is  not  quite  as 
obvious  when  the  decision  is  made  whether  or  not  to  allocate  K  * 
aircraft.  The  benefit  to  get  from  K-j*  to  K2*  is  equal  to  the  incremental 
capacity  reduction  multiplied  by  D. 

Benefit  =  62.23  X  7,500 

=  466,725. 
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Figure  5.  Example  Network  Capacity 
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This  benefit  is  compared  to  the  cost  of  allocating     K^*  -  K-|*    aircraft. 

Cost     =     20  X  30,000 

=     600,000. 

Thus  the  benefit  is  outweighed  by  the  cost.  Since  there  is:  na  allocation 
of  aircraft  greater  than  K2*  that  will  result  in  benefits  exceeding 
costs,  it  may  be  concluded  that  K*  =  41  represents  the  optima:!  number 
of  aircraft  to  assign  to  the  airstrike.  At  this  level  of  interdiction 
the  network  capacity  will  be  1056.06  tons  per  day.  This  is  a  reduction 
of  338.94  tons  per  day  with  a  resulting  benefit  of  Z  ,.542,050  dollars. 
The  cost  of  this  reduction  is  1,230,000  dollars.  The  cut  set:  that  will 
be  attacked  is  the  cut  set  corresponding  to  path  number  three  which 
contains  the  following  primal  arcs  :  (4,7);  (4,6);  (2,6);  (1,6);  and 
(1,3).  These  arcs  correspond  to  dynamic  programming  stages  1,2,3,4, 
and  5  respectively.  Looking  at  the  dynamic  programming  stages  the 
optimal  allocation  of  aircraft  is:  k^  7  =  9;  k«  6  =  6;  k~  6  =  7; 
k-^5  =  10;  and  k-|  3  =  9.  This  completes  the  solution. 
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IV.   DISCUSSION 

A.   PROPERTIES  OF  THE  SOLUTION  TECHNIQUE 

The  dynamic  programming  approach  taken  to  the  problem  guarantees 
that  the  solution  found  will  be  a  global  minimum  aver  the  feasible 
region.  The  integer  constraints  pose  no  problem.  En  fact,  the  integer 
restriction  limits  the  number  of  values  the  decision  variabiles:  may 
assume  and  allows  an  exact  solution  to  be  found.  Dynamic  programming 
also  provides  a  built  in  capability  for  sensitivity  analysis. 

The  convexity  of  the  damage  function  allowed  the  use  of  Fibonacci 
search  within  the  dynamic  program  resulting  in  a  tremendous  savings  in 
the  number  of  separate  calculations  made  in  each  dynamic  program.  With 
K  equal  to  100  the  reduction  was  on  the  order  of  10"'  times  the  number 
of  calculations  needed  for  exhaustive  search.  Larger  vaTues  of  K  will 
produce  savings  of  an  even  greater  magnitude.  The  execution  time 
required  for  the  sample  problem  was  15.06  seconds  on  an  IBM  360/67. 
Utilizing  Fibonacci  search  it  was  found  that  the  increase  in  execution 
time  for  larger  values  of  K  was  approximately  linear.  Execution  time 
was  also  roughly  linear  with  respect  to  the  total  number  of  dynamic 
programming  stages  required  (45  in  the  sample  problem).  From  the  above 
observations  the  amount  of  computer  time  required  for  larger  problems 
may  be  predicted.  For  example,  a  problem  with  15  cut  sets  in  S  averaging 
6  arcs  per  cut  set  would  require  90  dynamic  programming  stages.  If  200 
aircraft  were  available  a  reasonable  estimate  would  be  that  this  problem 
would  take  approximately  4  times  as  long  to  solve  as  the  sample  problem. 

A  further  reduction  in  the  number  of  calculations  required  may  be 
achieved  with  a  coarse  grid.  Aircraft  can  be  allocated  in  packages  of 
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five  and  the  constraining  cut  sets  determined:.     These  cut  sets  cair  then 
have  aircraft  reallocated  one  at  a  time  and  the  optimal   solution  found 
as  before.     This  approach  can  not  guarantee  that  the  correct  constraining 
cut  sets  will  be  selected,  but  if  they  are  the  solution  will  be  optimal. 

Dynamic  programming  allows  some  generalizations  to  be  made  in  the 
problem.     To  begin  with,  since  additive  stage  returns  aree  always: 
decomposable,  the  technique  places  no  restrictions  oar  the  damage^  functions. 
The  negative  exponential   damage  function  used  in  this  paper  has:  intuitive 
appeal   since  it  does  exhibit  diminishing  marginal   returns.     This  function 
also  contributes  to  computational  efficiency  since  its  convexity  allowed 
the  use  of  Fibonacci  search.     However,   if  actual  interdiction  data 
suggests  damage  functions  of  another  form,  the  problem  can  ssttIT  be- 
solved  with  somewhat  greater  expenditures  of  computer  time. 

Another  generalization  suggested  by  dynamic  programming  is:  to:  consider 
the  allocation  of  two  types  of  aircraft.     In  this  case  a  damage  function 
of  the  form 

mij(kij>hij)     =     ]ij  +wij  exP("bijkij   "  aijhij) 

might  be  assumed  with  k-jj,  1-jj,  and  w-j,-  defined  as  before,  h^  •   repre- 
senting the  number  of  aircraft  of  the  second  type  assigned  to  arc  (i,j), 
and  a^j  denoting  the  vulnerability  parameter  corresponding  to  the 
second  type  of  aircraft.     Dynamic  programming  may  again  be  used  to  solve 
the  problem,  but  two  state  and  two  decision  variables  are  required. 

Although  the  new  damage  function  preserves  convexity,  in  this  case 
the  series  of  minimizations  is  of  functions  of  two  variables  and 
Fibonacci   search  is  not  applicable.     A  minimization  problem  was  run  for 
a  hypothetical   cut  set  containing  five  arcs.     The  execution  time  required 
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for  solution  was  5.63  seconds  when  TO  aircraft  of  two  types:  were 

available;  with  19  aircraft  of  each  type,  the  time  required  was  32.79 

seconds;  and  when  25  of  each  type  aircraft  were  available,  over  a 

minute  of  computer  time  was  used.     To  deal  with  even  relatively  small 

networks  the  computer  time  requirements  would  become  prohibitive  if  it 

were  necessary  to  consider  larger  aircraft  availabilities.     To:  assign 

three  types  of  aircraft,  dynamic  programming  would:  require  three:  state 

and  three  decision  variables  and  the  technique  would  be  impractical 

even  for  small   problems. 

Another  application  of  the  dynamic  programming  approach  is:  in  a. 

modification  of  Nugent's  algorithm  [3].     This  modification  will   provide 

integer  solutions.     Nugent  presented  a  method  of  finding  non-integer 

allocations  of  resources  that  would  minimize  network  capacity  subject 

to        z  k.-_-  <  K      and     k . .  >  0   .     As  previously  discussed,  the 

U,J)eSi   1J  -  *T- 

objective  function  in  this  problem  is  convex  with  respect  to  k...  In 

'  j 

Nugent's  formulation  the  feasible  region  defined  by  the  constraints  is 
also  convex.     Therefore,  for  any  particular  cut  set  the  problem  is  a 
convex  non-linear  program  and  Kuhn-Tucker  theory  provides  conditions 
that  are  both  necessary  and  sufficient  for  a  global  minimum.     Nugent 
solves  these  Kuhn-Tucker  conditions  and  using  an  upper  bounding  technique 
arrives  at  the  cut  set  that  will  be  minimal  after  optimum  interdiction. 

In  the  modification  the  set  S  and  the  upper  bound  on  network  capacity 
are  found  as  before.     The  Kuhn-Tucker  conditions  are  then  solved  to  find 
non-integer  constrained  solutions  that  minimize  cut  set  capacity  for 
each  element  of  S.     The  minimum  of  these  solutions  represents  the  optimal 
so'lution  without  integer  constraints.     When  integer  constraints;  are  added 
this  minimum  represents  a  lower  bound  on  network  capacity.     The  cut  set 
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with  the  minimal   non-integer  solution  is  deleted  from  S  arret  becomes  the: 
subject  of  a  dynamic  program  to  find  an  integer  solution.     Lf  this 
integer  solution  is  less  than  the  non-integer  solutions  corresponding 
to  the  remaining  elements  of  S  it  is  optimal.     If  it  is  greater  than  some 
or  all  of  the  elements  of  S  it  represents  a  new,  smaller  upper  bound  on 
network,  capacity.     Any  elements  of  S  with  nan-integer  solutions  greater 
than  this  new  upper  bound  are  deleted  from  S.     From  the  remaining:  elements: 
of  S  the  cut  set  with  the  smallest  non-integer  solutforr  is  selected  from 
S.     Again  dynamic  programming  used  to  find  a  new  integer  solution.     The 
new  integer  solution  is  compared  with  the  old  integer  solution  and  the 
minimum  is  called  the  current  integer  solution.     The  current  integer 
solution  is  then  compared  to  the  remaining  non-integer  solutions  and 
the  process  is  repeated.     This  iterative  procedure  is  continued  until 
either  S  is  the  null   set  or  until   the  current  integer  solution  is  less 
than  or  equal   the  non-integer  solutions  corresponding  to  all   of  the 
remaining  elements  of  S.     In  either  case  the  current  integer  solution 
represents  the  optimal  solution  to  the  integer  constrained  problem. 

In  general,  if  the  number  of  aircraft  to  be  allocated  to  the  air- 
strike  is  known,  this  modification  is  more  efficient  that  using  dynamic 
programming  on  every  element  of  S  to  solve  the  minimization  problem. 
In  solving  the  example  problem  from  Nugent's  paper  it  was  necessary  to 
run  only  one  dynamic  program  and  with  exponential   damage  functions  the 
Kuhn-Tucker  conditions  are  easy  to  solve  relative  to  solving  a  dynamic 
program.     However,  this  modification  does  not  lend  itself  to  the 
sensitivity  analysis  on  K  that  is  necessary  when  the  number  of  aircraft 
to  be  assigned  to  the  strike  is  taken  to  be  a  decision  variable. 
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As  already  mentioned,  there  are  limitations  an  the:  technique 
presented.     One  difficulty  that  has  not  yet  been  discussed  is  in  the 
measurement  of  the  costs  and  benefits  of  aircraft  assignment.     In  this 
paper  the  problem  was  ignored  and  constant  dollar  values  of  C  and  D 
were  selected  arbitrarily.     This  problem  ts  important  .since:  the  selection 
of  C  and  D  determines  how  many  aircraft  wtTT  he:  assigned:  to  tte  strike. 
If  D  had  been  taken  to  be  10,000  dollars  per  torr  of  flow  reduced:  vice. 
7,500  and  the  rest  of  the  problem  remained:  unchanged,   the:  decision  would 
have  been  made  to  allocate  11  aircraft  in  a  strike  against  cut  set  one 
resulting  in  a  network  capacity  of  938.19  tons  per  day.     On  the:  other 
hand,  if  D  was  less  than  1,519  dollars  per  tan  of  flow  reduced:  the. 
solution  would  be  to  make  no  attack  against  the  network. 

B.      RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  possibility  of  deriving  damage  functions  from  actual   interdiction 
data  was  mentioned  earlier.     If  the  method  of  this  paper  were  to  be  put 
to  use  in  solving  a  real -world  interdiction  problem  some  verification 
of  the  damage  function  would  be  essential.     However,  due  to  the  sensi- 
tivity of  the  solution  to  both  costs  and  benefits,  the  measurement 
problem  associated  with  costs  and  benefits  should  receive  at  least  as 
much  attention  as  the  damage  function. 

Another  possibility  for  further  study  would  be  the  utilization  of 
the  model   described  in  this  paper  to  represent  real -world  problems  other 
than  aircraft  interdiction.     One  obvious  example  might  be  the  problem  of 
allocating  resources  to  the  improvement  of  a  highway  system.     In  this 
example  it  would  probably  be  relatively  easy  to  get  data  from  which 
to  derive  improvement  functions,  but  the  measurement  of  costs  and  benefits 
would  be  as  difficult  as  before. 
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The  model  presented  could  he  refined  try  assigning  different, values 
to  capacity  reduction  in  the  various  arcs  of  the  network.  The  objective 
then  would  be  to  minimize  the  maximum  value  of  flow  possible  in  the 
network  rather  than  to  minimize  network  capacity.  The  solution  technique 
presented  could  still  be  used.  A  further  refinement  might  be  to  consider 
not  only  arc  vulnerability,  but  also  the  repair  capability  of- the 
opponent.  This  would  require  capacity  reduction  to:  be^  taken  as  a 
function  of  time  as  well  as  aircraft  allocation  and:  would-  make  .the 
analysis  of  the  model  more  difficult.  Many  other  refinements  could  be 
made  in  order  to  make  the  model  more  representative  of  the:  real  world, 
but  in  general  the  increased  realtsm  gained:  would:  be  at.  the: expense  of 
increased  computational  effort. 
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V.      SUMMARY 


A  solution  procedure  has  been  developed  for  the  problem  of 
determining  the  optimal   allocation  of  aircraft  in  planning  an  airstrike 
against  a  transportation  network.     The  damage  function  for  arcs:  in  ther 
network  is  assumed  to  have  a  negative  exponential  form.     Tbimake^use- 
of  the  procedure  it  is  necessary  to  have  available  the  following-  infor- 
mation:    the  upper  and  lower  bounds  on  the  capacity  of  each  arc^  the 
vulnerability  parameter  for  each  arc,  the  number  of  aircraft  available 
for  assignment  to  the  airstrike,  the  cost  of  assigning  an  aircraft:  to 
the  strike,  and  the  benefit  resulting  from  network  capacity  reduction. 

In  the  solution  procedure  es/ery  cut  set  that  is  designated  a 
candidate  for  attack  is  the  subject  of  a  dynamic  program.     A  sensitivity 
analysis  is  performed  on  the  number  of  aircraft  to  be  assigned  and  this 
gives  the  network  capacity  after  optimal   interdiction  as  a  function  of 
the  number  of  aircraft  assigned  to  the  strike.     A  cost  benefit  analysis 
is  then  made  to  determine  the  largest  number  of  aircraft  that  can  be 
assigned  before  costs  of  further  allocation  begin  to  outweigh  the 
benefits  resulting  from  that  allocation. 

At  the  end  of  the  procedure  the  solution  consists  of  the  following: 
the  number  of  aircraft  to  assign  to  the  airstrike,  the  cut  set  that  will 
be  attacked,  the  number  of  aircraft  to  allocate  to  each  arc  of  the  cut 
set  chosen,  and  the  capacity  of  the  network  after  this  assignment  of 
aircraft. 
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COMPUTER  OUTPUT 
DP    SOLNS    CUT    SET    1 


ACFT  AVAIL 

CAPACITY 

0 

1500.0000 

1 

1472.2615 

2 

1450.4417 

3 

1430.5337 

4 

1413.3696 

5 

1397. 8652 

6 

1384.3635 

7 

1372.2888 

8 

1360.5256 

9 

1349.2236 

10 

1338.3647 

11 

1327.7441 

12 

1317.3113 

13 

1307.2874 

14 

1297.5332 

15 

1267.9023 

16 

1278.4985 

17 

1269.2200 

18 

1259.9668 

19 

1251.0762 

20 

1242.2502 

21 

1233. 70E5 

22 

1225.3132 

23 

1216.9587 

24 

1208.7517 

25 

1200.7656 

26 

11S2.88C6 

27 

1165.2842 

28 

1177.7083 

29 

1170.3845 

30 

1163.1057 

31 

1155.8796 

32 

1148.8862 

33 

1142.0127 

34 

1135.2935 

35 

1128.7214 

36 

1122.1831 

37 

1115.7273 

38 

1109.5078 

39 

1103.3052 

40 

1C<57.3457 

41 

1091.4297 

42 

1085.7039 

43 

1C80.0000 

44 

1074.3726 

45 

1068.8713 

46 

1063.5181 

47 

1058.2324 

48 

1053.0627 

49 

1047.9707 

CHANGE    IN    CAPACITY 


0.0 

-27.7385 

-21.8198 

-19.9080 

-17.1641 

-15.5044 

-13.5017 

-12.0747 

-11.7632 

-11.3020 

-10.8589 

-10.6206 

-10.4329 

-10.0239 

-9.7542 

-9.6309 

-9.4038 

-9.2786 

-9.2532 

-8.8906 

-8.8259 

-8.5417 

-8.3953 

-8.3545 

-8.2070 

-7.9861 

-7.8850 

-7.5964 

-7.5759 

-7.3237 

-7.2788 

-7.2261 

-6.9934 

-6.8735 

-6.7192 

-6.5720 

-6.5383 

-6.4558 

-6.2195 

-6.2026 

-5.9595 

-5.9160 

-5.7258 

-5.7039 

-5.6274 

-5.5012 

-5.3533 

-5.2856 

-5.1697 

-5.0920 
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ACFT   AVAIL 

CAPACITY 

50 

1042,8923 

51 

1038.0132 

52 

1033.1694 

53 

1028.4817 

54 

1023.8743 

55 

1019.3701 

56 

1014.9280 

57 

1010.5452 

58 

1006.2178 

59 

1002.0488 

60 

997.8911 

61 

993.8245 

62 

989.8296 

63 

985.8640 

64 

982.0259 

65 

978.2537 

66 

974.5659 

67 

970.9778 

68 

967.4346 

69 

963.9751 

70 

960.5618 

71 

957.1577 

72 

953.8872 

73 

950.6404 

74 

947.4414 

75 

944.2991 

76 

941.2104 

77 

938.1914 

78 

935.2534 

79 

932.3528 

80 

929.5581 

81 

926.7710 

82 

924.0769 

83 

921.3992 

84 

918.7407 

85 

916.1680 

86 

913.6394 

87 

911.1230 

88 

908.6511 

89 

906.2458 

90 

903.8708 

91 

901.5830 

92 

899.3010 

93 

897.1086 

94 

894.9324 

95 

£92.8259 

96 

890.7275 

97 

888.6572 

98 

886.6335 

99 

884.6541 

100 

882.6848 

CHANGE    IN    CAPACITY 


-5.0784 
-4.8792 
-4.8435 
-4.6879 
-4.6075 
-4.5040 
-4.4421 
-4.3828 
-4.3274 
-4.1691 
-4.1577 
-4.0666 
-3.9947 
-3.9657 
-3.8381 
-3.7723 
-3.6876 
-3.5883 
-3.5430 
-3.4595 
-3.4133 
-3.4041 
-3.270  6 
-3.2469 
-3.1989 
-3.1423 
-3.0885 
-3.0191 
-2.9  3  79 
-2.9008 
-2.7946 
-2.7870 
-2.6942 
-2.6777 
-2.6583 
-2.5727 
-2.5287 
-2.5164 
-2.4719 
-2.4053 
-2.3749 
-2.2880 
-2.2818 
-2.1924 
-2.1765 
-2.1064 
-2.0983 
-2.0703 
-2.0238 
-1.9794 
-1.9693 
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STAGE    NUMBER 


1ST ATE  IDEC      I  STATE  I  DEC      I  STATE  IDEC 


0 

0 

1 

1 

2 

2 

3 

3 

4 

4 

5 

5 

6 

6 

7 

7 

a 

8 

9 

9 

10 

10 

ii 

LI 

12 

12 

13 

13 

14 

14 

15 

15 

16 

16 

17 

17 

18 

18 

19 

19 

20 

20 

21 

21 

22 

22 

23 

23 

24 

24 

25 

25 

26 

26 

27 

27 

28 

28 

29 

29 

30 

30 

31 

31 

32 

32 

33 

33 

34 

34 

35 

35 

36 

36 

37 

37 

38 

38 

39 

39 

40 

40 

41 

41 

42 

42 

43 

43 

44 

44 

45 

45 

46 

46 

47 

47 

48 

48 

49 

49 

50 

50 

51 

51 

52 

52 

53 

53 

54 

54 

55 

55 

56 

56 

57 

57 

58 

58 

59 

59 

60 

60 

61 

61 

62 

62 

63 

63 

64 

64 

65 

65 

66 

66 

67 

67 

68 

68 

69 

69 

70 

70 

71 

71 

72 

72 

73 

73 

74 

74 

75 

75 

76 

76 

77 

77 

78 

78 

79 

79 

80 

80 

81 

81 

82 

82 

83 

83 

84 

84 

85 

85 

86 

86 

87 

87 

88 

88 

89 

89 

90 

90 

91 

91 

92 

92 

93 

93 

94 

94 

95 

95 

96 

96 

97 

97 

98 

98 

99 

99 

100 

100 
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STAGE    NUMBER 


I  STATE  IDEC      I  STATE  IDEC      ISTATE  I.DEC 


0 

0 

1 

0 

2 

a 

3 

0 

4 

0 

5 

L 

6 

2 

7 

3 

a 

33 

9 

4 

10 

5 

n 

6c 

12 

7 

13 

8 

14 

* 

15 

9 

16 

10 

17 

111. 

18 

12 

19 

13 

20 

14 

21 

15 

22 

15 

23 

16 

24 

17 

25 

18 

26 

19 

27 

20 

28 

21 

29 

21 

30 

22 

31 

23 

32 

24 

33 

25 

34 

26 

35 

27 

36 

27 

37 

28 

38 

29 

39 

30 

40 

31 

41 

32 

42 

33 

43 

33 

44- 

34- 

45 

35 

46 

36 

4-7 

3T 

48 

38 

49 

39 

50 

39^ 

51 

40 

52 

41 

53 

42 

54 

43 

55 

44 

56 

45 

57 

45 

58 

46 

59 

47 

60 

48 

61 

49 

62 

50 

63 

51 

64 

51 

65 

52 

66 

53 

67 

54 

68 

55 

69 

56 

70 

57 

71 

57 

72 

58 

73 

59 

74 

60 

75 

61 

76 

62 

77 

63 

78 

63 

79 

64 

80 

65 

81 

66 

82 

67 

83 

68 

84 

69 

85 

69 

86 

70 

87 

71 

88 

72 

89 

73 

90 

74 

91 

75 

92 

75 

93 

76 

94 

77 

95 

78 

96 

79 

97 

80 

98 

81 

99 

81 

100 

82 
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STAGE    NUMBER      3 


1ST ATE  IDEC      ISTATE  IDEC      I  STATE  IDEC 


0 

0 

1 

0 

2 

0 

3 

1 

4 

1 

5 

2 

6 

2 

7 

3 

8 

3? 

9 

3 

10 

3 

11 

3 

12 

3 

13 

3 

14 

3 

15 

4 

16 

4 

17 

4 

18 

4 

19 

4 

20 

4 

21 

4 

22 

4 

23 

5 

24 

5 

25 

5 

26 

5 

27 

5 

28 

5 

29 

5 

30 

5 

31 

5 

32 

6 

33 

6 

34 

6 

35 

6r 

36 

6 

37 

6 

38 

6 

39 

6 

40 

7 

41 

7 

42 

7 

43 

7 

44 

7 

45 

7 

46 

7 

47 

7 

48 

8 

49 

8 

50 

a 

51 

8 

52 

8 

53 

8 

54 

8 

55 

8 

56 

9 

57 

9 

58 

9 

59 

9 

60 

9 

61 

9 

62 

9 

63 

9 

64 

9 

65 

10 

66 

10 

67 

10 

68 

10 

69 

10 

70 

10 

71 

10 

72 

10 

73 

11 

74 

LI 

75 

11 

76 

11 

77 

11 

78 

11 

79 

11 

80 

11 

81 

12 

82 

12 

83 

12 

84 

12 

85 

12 

86 

12 

87 

12 

88 

12 

89 

13 

90 

13 

91 

13 

92 

13 

93 

13 

94 

13 

95 

13 

96 

13 

97 

13 

98 

14 

99 

14 

100 

14 
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STAGE    NUMBER 


I  STATE  IDEC      I  STATE  IDEC      1ST  ATE  IDEC 


0 

0 

1 

0 

2 

0 

3 

0 

4 

0 

5 

0 

6 

0 

7 

0 

8 

0 

9 

0 

10 

0 

11 

0 

12 

0 

13 

0 

14 

L 

15 

1 

16 

1 

17 

2 

18 

2 

19 

2 

20 

3 

21 

3 

22 

4 

23 

4 

24 

4 

25 

5 

26 

5 

27 

6 

28 

6 

29 

6 

30 

6 

31 

7 

32 

7 

33 

8 

34 

8 

35 

8 

36 

9 

37 

9 

38 

10 

39 

10 

40 

10 

41 

LI 

42 

11 

43 

11 

44 

12 

45 

12 

46 

13 

47 

13 

48 

13 

49 

14 

50 

14 

51 

14 

52 

15 

53 

15 

54 

16 

55 

16 

56 

16 

57 

17 

58 

17 

59 

18 

60 

18 

61 

18 

62 

18 

63 

19 

64 

19 

65 

20 

66 

20 

67 

21 

68 

21 

69 

21 

70 

22 

71 

22 

72 

22 

73 

23 

74 

23 

75 

23 

76 

24 

77 

24 

78 

25 

79 

25 

80 

26 

81 

26 

82 

26 

83 

26 

84 

27 

85 

27 

86 

28 

87 

28 

88 

28 

89 

29 

90 

29 

91 

30 

92 

30 

93 

30 

94 

31 

95 

31 

96 

31 

97 

32 

98 

32 

99 

32 

100 

33 

39 


DP    SOLNS   CUT    SET    3 


AVAIL 

CAPACITY 

0 

1535.0000 

1 

1497.0632 

2 

1469.3250 

3 

1443.3813 

4 

1421.5618 

5 

1401.7942 

6 

1384.0525 

7 

1366.6409 

8 

1349.4768 

9 

1332.9653 

10 

1317.1670 

11 

1302.1807 

12 

1288.3894 

13 

1274.8879 

14 

1261.9890 

15 

1249.3105 

16 

1237.1775 

17 

1225.6582 

18 

1214.5559 

19 

1203.9353 

20 

1193.7605 

21 

1184.1384 

22 

1174.5828 

23 

1166.2283 

24 

1157.9309 

25 

1149.706  3 

26 

1141.5410 

27 

1133.5042 

28 

1126.4250 

29 

1119.7122 

30 

1113.1401 

31 

1106.5872 

32 

1100.4941 

33 

1094.8201 

34 

1089.2129 

35 

1083.9541 

36 

1078.7097 

37 

1073.5400 

38 

1068.8567 

39 

1064.3428 

40 

1060.1226 

41 

1056.0559 

42 

1052.1440 

43 

1048.2588 

44 

1044.3787 

45 

1040.9917 

46 

1037.6477 

47 

1034.3801 

48 

1031.1814 

49 

1028.3032 

CHANGE    IN    CAPACITY 


0.0 

-37.9368 

-27.738  3 

-25.9436 

-21.8196 

-19.7676 

-17.7417 

-17.4116 

-17.1641 

-16.5115 

-15.7983 

-14.9863 

-13.7913 

-13.5015 

-12.8989 

-12.6785 

-12.1331 

-11.5193 

-11.1021 

-10.6208 

-10.1748 

-9.6219 

-9.5556 

-8.3547 

-8.2972 

-8.2246 

-8.1655 

-8.0369 

-7.0790 

-6.7130 

-6.5720 

-6.5529 

-6.0930 

-5.6741 

-5.6072 

-5.2589 

-5.2442 

-5.1697 

-4.6835 

-4.5138 

-4.2203 

-4.0666 

-3.9120 

-3.8850 

-3.8803 

-3.3869 

-3.3439 

-3.2675 

-3.1989 

-2.8781 
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ACFT  AVAIL 

CAPACITY 

50 

1025.5740 

51 

1022.8560 

52 

1020.2021 

53 

1017.6858 

54 

1015.2087 

55 

1012.9290 

56 

1010.7476 

57 

1CC8.6155 

58 

1006.6360 

59 

1004.7319 

60 

1002.8967 

61 

1001.0820 

62 

999.3315 

63 

997.7410 

64 

996.1616 

65 

994.6045 

66 

993.1997 

67 

991.8401 

68 

990.5115 

69 

989.2705 

70 

988.0457 

71 

986.8755 

72 

985.7483 

73 

984.6384 

74 

983.6313 

75 

982.6680 

76 

981.7410 

77 

980.8362 

78 

979.9695 

79 

979.1208 

80 

978.3464 

81 

977.5886 

82 

976.8425 

83 

976.1165 

84 

975.4697 

85 

974.8276 

86 

974.2314 

87 

973.6487 

88 

973.0684 

89 

972.5156 

90 

971. 9753 

91 

971.4998 

92 

971.0308 

93 

970.5630 

94 

970.1121 

95 

969.7024 

96 

969.3057 

97 

968.9287 

98 

968.5535 

99 

968.1846 

100 

967.8320 

CHANGE    IN    CAPACITY 


-2.7293 
-2.7181 
-2.6536 
-2.5164 
-2.4772 
-2.2797 
-2.1813 
-2.1321 
-1.9794 
-1.9042 
-1.8352 
-1.8147 
-1.7505 
-1.5905 
-1.5795 
-1.5571 
-1.4048 
-1.3595 
-1.3285 
-1.2410 
-1.2249 
-1.1702 
-1.1274 
-1.1096 
-1.0072 
-0.9635 
-0.9269 
-0.9048 
-0.8669 
-0.8487 
-0.7742 
-0.7579 
-0.7461 
-0.7261 
-0.6466 
-0.6422 
-0.5962 
-0.5827 
-0.5804 
-0.5527 
-0.5401 
-0.4757 
-0.4690 
-0.4676 
-0.4511 
-0.4095 
-0.3969 
-0.3768 
-0.3753 
-0.3689 
-0.3524 
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STAGE    NUMBER 


ISTATE  IDEC      ISTATE  IDEC       I  STATE  IDEC 


0  0  112  2 

3  3  4  4  5  5 

6  6  7  T  8:  8 

9  g  io  ia  n:  n 

12  12  13  13  14  14 

15  15  16  16  IT  17 

18  18  19  19  20  20 

21  21  22  22  23  23 

24  24  25  25  26  26 

27  27  28  28  29  29 

30  30  31  31  32  32 

33  33  34  34  3.5  35 

36  36  37  37  38  38 

39  39  40  40  4-11  41 

42  42  43  43  44-  44 

45  45  46  46  4T  47 

48  48  49  49  50  50 

51  51  52  52  53  53 

54  54  55  55  56  56 

57  57  5  8  5  8  59  59 

60  60  61  61  62  62 

63  63  64  64  65  65 

66  66  67  67  68  68 

69  69  70  70  71  71 

72  72  73  73  74  74 

75  75  76  76  77  77 

78  78  79  79  80  80 

81  81  82  82  83  83 

84  84  85  85  86  86 

87  87  88  88  89  89 

90  90  91  91  92  92 

93  93  94   '     94  95  95 

96  96  97  97  98  98 

99  99  100  100 
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STAGE    NUMBER 


I  STATE  IDEC      ISTATE  IDEC      1ST  ATE  IDEC 


0 

0 

1 

1 

2 

1 

3 

2 

4 

2 

5 

3 

6 

3 

7 

3 

a 

4- 

9 

4 

10 

4 

n 

5- 

12 

5 

13 

6 

14 

6 

15 

6 

16 

7 

17 

7" 

18 

8 

19 

8 

20 

8 

21 

9 

22 

9 

23 

10 

24 

10 

25 

10 

26 

11 

27 

11 

28 

11 

29 

12 

30 

12 

31 

13 

32 

13 

33 

13 

34 

14 

35 

1.4 

36 

15 

37 

15 

38 

15 

39 

16 

40 

16 

41 

L6 

42 

17 

43 

17 

4A 

1.8 

45 

18 

46 

18 

47 

1.9 

48 

19 

49 

20 

5a 

2a 

51 

20 

52 

21 

53 

21 

54 

22 

55 

22 

56 

22 

57 

23 

58 

23 

59 

23 

60 

24 

61 

24 

62 

25 

63 

25 

64 

25 

65 

26 

66 

26 

67 

27 

68 

27 

69 

27 

70 

28 

71 

28 

72 

28 

73 

29 

74 

29 

75 

30 

76 

30 

77 

30 

78 

31 

79 

31 

80 

32 

81 

32 

82 

32 

83 

33 

84 

33 

85 

34 

86 

34 

87 

34 

88 

35 

89 

35 

90 

35 

91 

36 

92 

36 

93 

37 

94 

37 

95 

37 

96 

38 

97 

38 

98 

39 

99 

39 

100 

39 
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STAGE  NUMBER   3 


ISTATE  IDEC      ISTATE  IDEC      I  STATE  IDEC 


0 

0 

1 

0 

2 

0 

3 

0 

4 

0 

5 

0 

6 

0 

7 

1 

8 

1 

9 

2 

10 

2 

11 

2 

12 

3 

13 

3 

14 

3 

15 

4 

16 

4 

17 

5 

18 

5 

19 

6 

20 

6 

21 

7 

22 

7 

23 

7 

24 

8 

25 

8 

26 

9 

27 

9 

28 

9 

29 

10 

30 

10 

31 

10 

32 

11 

33 

11 

34 

12 

35 

12 

36 

12 

37 

13 

38 

13 

39 

14 

40 

14 

41 

14 

42 

15 

43 

15 

44 

16 

45 

16 

46 

16 

47 

17 

48 

17 

49 

18 

50 

18 

51 

19 

52 

1<3 

53 

19 

54 

20 

55 

20 

56 

21 

57 

21 

58 

21 

59 

22 

60 

22 

61 

22 

62 

23 

63 

23 

64 

24 

65 

24 

66 

24 

67 

25 

68 

25 

69 

26 

70 

26 

71 

26 

72 

27 

73 

27 

74 

28 

75 

28 

76 

28 

77 

29 

78 

29 

79 

30 

80 

30 

81 

31 

82 

31 

83 

31 

84 

32 

85 

32 

86 

33 

87 

33 

88 

33 

89 

34 

90 

34 

91 

35 

92 

35 

93 

35 

94 

36 

95 

36 

96 

36 

97 

37 

98 

37 

99 

38 

100 

38 
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STAGE  NUMBER 


I  STATE     IDEC   I  STATE     IDEC   1ST  ATE     IDEC 


0  0  10  2  0 

3  0  4  0  5  0 

6  17  18  1 

9  2  10  2  11  3 

12  3  13  3  14  4 

15  4  16  4  17  5 

18  5  19  5  20  6 

21  6  22  7  23  7 

24  7  25  8  26  8 

27  8  28  9  29  9 

30  10  31  10  32  10 

23  11  34  11  35  11 

36  12  37  12  38  13 

39  13  40  13  41  13 

42  14  43  14  44  15 

45  15  46  16  47  16 

48  16  49  17  50  17 

51  17  52  18  53  18 

54  18  55  19  56  19 

57  20  58  20  59  20 

60  21  61  21  62  21 

63  22  64  22  65  23 

66  23  67  23  68  23 

69  24  70  25  71  25 

72  25  73  26  74  26 

75  26  76  26  77  27 

78  28  79  28  80  28 

81  28  82  29  83  29 

84  29  85  30  86  30 

87  31  88  31  89  31 

90  32  91  32  92  33 

93  33  94  33  95  33 

96  34  97  34  98  35 

99  35  100  36 
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STAGE  NUMBER   5 


ISTATE 


IDEC   ISTATE 


IDEC   ISTATE 


IDEC 


0 

3 
6 
9 

12 

15 

18 

21 

24 

27 

30 

33 

36 

39 

42 

45 

48 
51 
54 
57 

60 

63 

66 

69 

72 

75 

78 

81 

84 

87 

90 

93 

96 

99 


0 
0 

1 
2 
3 
3 
4 
5 
5 
6 
7 
7 
8 
9 
10 
10 
11 
12 
12 
13 
14 
15 
15 
16 
16 
17 
18 
19 
20 
20 
21 
21 
22 
23 


1 
4 
7 

10 

13 

16 

19 

22 

25 

28 

31 

34 

37 

40 

43 

46 

49 

52 

55 

58 

61 

64 

67 

70 

73 

76 

79 

82 

85 

88 

91 

94 

97 

100 


0 
0 

1 
2 
3 
3 
4 
5 
5 
6 
7 
8 
8 
9 
10 
10 
11 
12 
13 
13 
14 
15 
15 
16 
17 
18 
18 
19 
20 
20 
21 
22 
23 
23 


2 

5 

a 
ii 

14 

17 

20 

23 

26 

29 

32 

35 

38 

41 

44 

47 

50 

53 

56 

59 

62 

65 

68 

71 

74 

77 

80 

83 

86 

89 

92 

95 

98 


Q 

I 
I 
2 
3 
4 
4 
5 
5 
7 
7 
8 
9 
9 
10 
II 
12 
12 
13 
14 
14 
15 
16 
16 
17 
18 
19 
19 
20 
20 
21 
22 
23 
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DP    SOLNS   CUT    SET    4 


ACFT  AVAIL 

CAPACITY 

0 

1395.0000 

1 

1375.2324 

2 

1357.8208 

3 

1341.3096 

4 

1326.3232 

5 

1312.5320 

6 

1298.9368 

7 

1286.0381 

8 

1274.5186 

9 

1263.3877 

10 

1252.2856 

11 

1242.6638 

12 

1233.1082 

13 

1223.9949 

14 

1215.7703 

15 

1207.7334 

16 

1200.2722 

17 

1193.1934 

18 

1186.480  2 

19 

1160.3716 

20 

1174.2786 

21 

1168.6714 

22 

1163.4272 

23 

1158.4258 

24 

1153.7424 

25 

1149.2285 

26 

1145.1338 

27 

1141.2219 

28 

1137.3367 

29 

1133.9841 

30 

1130.6404 

31 

1127.3728 

32 

1124.4946 

33 

1121.7500 

34 

1119.0205 

35 

1116.5435 

36 

1114.2637 

37 

1112.0164 

38 

1109.8843 

39 

11C7.9800 

40 

1106.1401 

41 

1104.3049 

42 

1102.7146 

43 

1101.1350 

44 

1C99.6287 

45 

1C98.2690 

46 

1096.9407 

47 

1C95.7073 

48 

1094.5371 

49 

1093.4275 

CHANGE    IN   CAPACITY 


0.0 

-19.7676 
-17.4116 
-16.5112 
-14.9861 
-13.7914 
-13.5952 
-12.8988 
-11.5195 
-11.1308 
-11.1021 
-9.6219 
-9.5556 
-9.1131 
-8.2246 
-8.0369 
-7.4612 
-7.0790 
-6.7130 
-6.1087 
-6.0930 
-5.6072 
-5.2442 
-5.0014 
-4.6835 
-4.5138 
-4.0948 
-3.9120 
-3.8850 
-3.3525 
-3.3439 
-3.2675 
-2.8781 
-2.7448 
-2.7293 
-2.4772 
-2.2797 
-2.2473 
-2.1322 
-1.9042 
-1.8399 
-1.8352 
-1.5905 
-1.5795 
-1.5064 
-1.3595 
-1.3285 
-1.2333 
-1.1702 
-1.1096 
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T    AVAIL 

CAPACITY 

50 

1092.4177 

51 

1091.4106 

52 

1090.4836 

53 

1C89.6169 

54 

1088.7900 

55 

1C88.0159 

56 

1087.2698 

57 

1086.5930 

58 

1065.9463 

59 

1085.3042 

60 

1084.7500 

61 

1C64.1973 

62 

1083.6570 

63 

1083.1814 

64 

1C82.7275 

65 

1082.2764 

66 

1C81.8669 

67 

1081.4902 

68 

1081.1187 

69 

1C80.7661 

70 

1080.4514 

71 

1080.1472 

72 

1C79.8440 

73 

1079.5811 

74 

1079. 32C1 

75 

1079.0710 

76 

1078.8462 

77 

1078.6267 

78 

1078.4229 

79 

1C78.2292 

80 

1C78.0459 

81 

1C77.8789 

82 

1077.7126 

83 

1C77.5593 

84 

1077.4160 

85 

1077.2793 

86 

1C77.1514 

87 

1077. C281 

88 

1076.9163 

89 

1076.8093 

90 

1076.7031 

91 

1C76.6116 

92 

1076.5203 

93 

1076.4309 

94 

1076.3523 

95 

1076.2773 

96 

1076.2026 

97 

1C76.1350 

98 

1076.0728 

99 

1C76.0U2 

100 

1C75.9531 

CHANGE    IN    CAPACITY 


-1.0098 
-1.0072 
-0.9269 
-0-8669 
-0-8267 
-0-7742 
-0-7461 
-0 -6769 
-0.6466 
-0.-6422 
-0-5542 
-0.5527 
-0.5401 
-0»4757 
-0.4537 
-0.4511 
-0.409  5 
-0-376  8 
-0-3715 
-0-3  524 
-0-3148 
-0-3041 
-0-30  3  3 
-0-2629 
-0.2611 
-0.2490 
-0-2247 
-0-2196 
-0.2039 
-0.1934 
-0.1834 
-0.1669 
-0.1665 
-0.1532 
-0.1433 
-0.1367 
-0.1280 
-0.1233 
-0.1119 
-0.1069 
-0.1062 
-0.0916 
-0.0914 
-0.0893 
-0.0786 
-0.0750 
-0.0746 
-0.0677 
-0.0623 
-0.0614 
-0.0583 
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STAGE  NUMBER 


ISTATE  IDEC      ISTATE  IOEC      ISTATE  IDEC 


0 

0 

1 

1 

2 

Z 

3 

3 

4 

4 

5 

5 

6 

6 

7 

7 

8 

8: 

9 

9 

10 

10 

11 

11L 

12 

12 

13 

13 

14 

14 

15 

15 

16 

16 

17 

IT 

18 

18 

19 

19 

20 

ZC 

21 

21 

22 

22 

23 

23 

24 

24 

25 

25 

26 

26 

27 

27 

28 

28 

29 

29 

30 

30 

31 

31 

32 

32 

33 

33 

34 

34 

35 

35 

36 

36 

37 

37 

38 

38 

39 

39 

40 

40 

41 

4L 

42 

42 

43 

43 

4-4 

44- 

45 

45 

46 

46 

47 

4T 

48 

48 

49 

49 

50 

50 

51 

51 

52 

52 

53 

53 

54 

54 

55 

55 

56 

56 

57 

57 

58 

58 

59 

59 

60 

60 

61 

61 

62 

62 

63 

63 

64 

64 

65 

65 

66 

66 

67 

67 

68 

68 

69 

69 

70 

70 

71 

71 

72 

72 

73 

73 

74 

74 

75 

75 

76 

76 

77 

77 

78 

78 

79 

79 

80 

80 

81 

81 

82 

82 

83 

83 

84 

84 

85 

85 

86 

86 

87 

87 

88 

88 

89 

89 

90 

90 

91 

91 

92 

92 

93 

93 

94 

94 

95 

95 

96 

S6 

97 

97 

98 

98 

99 

99 

100 

100 
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STAGE    NUMBER 


ISTATE  IDEC      ISTATE  IDEC      ISTATE  IDEC 


0  0  1  1  2  2 

3  2  4  3  5  3. 

6  4  7  5  a  5-: 

9  6  10  6  II  T 

12  7  13  8  14  9- 

15  9  16  10  17  10 

18  11  19  II  20  12 

21  13  22  13  23  14 

24  14  25  15  26  15 

27  16  28  17  29  IT 

30  18  31  18  32  19 

33  19  34  20  35  21 

36  21  37  22  38  22 

39  23  40  23  4-1  24 

42  25  43  25  44  26 

45  26  46  2  7  47  ZT 

48  28  49  29  50  Z9 

51  30  52  3  0  53  31. 

54  31  55  32  56  33 

57  33  58  34  59  34 

60  35  61  35  62  36 

63  27  64  37  65  38 

66  38  67  39  68  39 

69  40  70  41  71  41 

72  42  73  42  74  43 

75  43  7  6  44  77  45 

78  45  79  46  80  46 

81  47  82  47  83  48 

84  49  85  49  86  50 

87  50  88  51  89  51 

90  52  91  53  92  53 

93  54  94  54  95  55 

96  55  97  56  98  57 

99  57  100  58 
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STAGE    NUMBER      3 


STATE 

IDEC 

I  STATE 

IDEC   I 

STATE 

IDEC 

0 

0 

1 

1 

2 

1 

3 

2 

4 

2 

5 

3 

6 

3 

7 

3 

8 

4 

9 

4 

10 

4 

11 

5 

12 

5 

13 

5 

14 

5 

15 

6 

16 

6 

17 

6 

18 

7 

19 

7 

20 

T 

21 

8 

22 

8 

23 

8 

24 

9 

25 

9 

26 

9 

27 

10 

28 

10 

29 

10 

30 

10 

31 

11 

32 

LI 

33 

11 

34 

12 

35 

LZ 

36 

13 

37 

13 

38 

L3 

39 

14 

40 

14 

4-1 

L4 

42 

15 

43 

15 

44 

L5 

45 

15 

46 

16 

47 

L6l 

48 

16 

49 

17 

50 

LT 

51 

17 

52 

18 

53 

18 

54 

18 

55 

19 

56 

19 

57 

19 

58 

20 

59 

20 

60 

20 

61 

20 

62 

21 

63 

21 

64 

21 

65 

22 

66 

22 

67 

23 

68 

23 

69 

23 

70 

24 

71 

24 

72 

24 

73 

25 

74 

25 

75 

25 

76 

25 

77 

26 

78 

26 

79 

26 

80 

27 

81 

27 

82 

27 

83 

28 

84 

28 

85 

28 

86 

29 

87 

29 

88 

29 

89 

ao 

90 

30 

91 

30 

92 

30 

93 

31 

94 

31 

95 

31 

96 

32 

97 

32 

98 

33 

99 

33 

100 

33 
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C  COMPUTER    PROGRAM 

C 

c 

C  THIS    PROGRAM    IS    DESIGNED   TO    FIND   THE    MINIMUM    OF    A 

C  SEQUENCE    OF    FUNCTIONS.       EACH    FUNCTION    IN    THE    SEQUENCE 

C  IS    A    SUM    OF    NEGATIVE    EXPONENTIAL    FUNCTIONS    OF    THE 

C  FCRM    BLO+W*EXP(-B*IDEC).       IDEC    REPRESENTS    THE    AMOUNT 

C  OF    RESOURCE    ALLOCATED    TO    REDUCE    THE    VALUE    CF    A 

C  PARTICULAR    NEGATIVE    EXPONENTIAL    AND    IS    THE    DECISION 

C  VARIABLE.       THE     MINIMIZATION    IS    SUBJECT    TO    CONSTRAINT 

C  THAT    THE    SUM    OF    THE    IDEC'S    FOR    EACH    FUNCTION    IN    THE 

C  SEQUENCE    IS    LESS    THAN    OR    EQUAL    TO    K    WHERE    K 

C  REPRESENTS    THE     TOTAL    RESOURCE    AVAILABILTY.       DYNAMIC 

C  PROGRAMMING    IS    THE    SOLUTION    TECHNIQUE    USED. 

C 

C  THE    INPUT    PARAMETERS    FOR    THIS    PROGRAM    ARE    AS   FCELLOWS"  - 

C 

C  NCUT-NUMBER    OF    FUNCTIONS    TO    BE    MINIMIZED.. 

C  K    -    MAXIMUM   AMOUNT    OF    RESOURCE    AVAILABLE. 

C  UPBND    -    PREDETERMINED    UPPER    BOUND. 

C  N   -    NUMBER    OF    EXPONENTIALS    IN    PARTICULAR    FUNCTION. 

C  B,W,BLO    -    CONSTANTS    ASSOC    WITH    EACH    EXPONENTIAL. 

C 

C  THE  OUTPUT  FROM  THIS  PROGRAM  IS  AS  FOLLOWS  - 

C  (FOR  EACH  FUNCTION  N  DYNAMIC  PROGRAMMING  STAGES 

C  ARE  REQUIRED.   FOR  EACH  STAGE  THE 

C  OPTIMAL  VALUE  OF  THE  DECISION  VARIABLE  IS  PRINTED 

C  FOR  EACH  POSSIBLE  VALUE  OF  THE  STATE  VARIABLE..) 

C 

C  I  STATE  -  VALUE  OF  STATE  VARIABLE. 

C  IDEC  -  OPTIMAL  VALUE  FOR  DECISION  VARIABLE. 

C 

C  (AT  THE  END  OF  THE  STAGES  THE  SOLUTION  IS 

C  PRINTED  FOR  EACH  FUNCTION) 

C         AIRCRAFT  AVAIL  -  AMOUNT  OF  RESOURCE  AVAILABLE. 
C         CAPACITY  -  SOLUTION  TO  MINIMIZATION  PROBLEM. 
C         DELCAP  -  INCREMENTAL  IMPROVEMENT  IN  SOLUTION. 
C 

DIMENSION  FCN( 200) .FCNNEWC 200) ,CAP(200) 

DIMENSION  I  PLANE (200) ,KPL AN E( 200 ) T DELCAP ( 200) 

DIMENSION    IFIBNQC20) 

DATA    I FIBNO/ 1,2.4,7,1 2, 20,33,54, 8 8, 143 ,232,376, 609, 
1986,15  96,2  5  83,4180,6764,10945,17710/ 

C  READ    #    PATHS     IN    NETWORK,     #    PLANES    AVAIL 

C  &    LUB    ON    CAPACITY 

C 

READ(5,40)    NCUT,K, UPBND 
40    F0RMAT(2(I  3)  ,F9.4) 

DO    1200     IPATH=1,NCUT 

DATA    FCN/20C*0.0/,FCNNEW/200*0.0/,CAP/200*0.0/ 

DATA    IPLANE/200*0/,KPLANE/200*0/ 

DATA    DELCAP/200*0.0/ 
C 
C  READ    #   ARCS    IN    CUT    SET 


C 


C 


READ(5,50)    N 
50    F0RMAT(I2) 
M  =  K  +  1 


C 

C  READ    IN    ARC    PARAMETERS    (1ST    ARC) 

C 

READ(5,60)     B,W,BLO 
60    F0RMAT(3(F10.6) ) 

BLOSUM=BLO 
C 
C  DP    STAGE    1 


DO    100    I=1,M 
ARG=-B*(I-1) 
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FCN(I) =W*EXP(ARG) 
IPLANE(  I  )=  I- 1 
KPLANEU  )=IPLANE(I  ) 
100    CONTINUE 
ISTAGE=1 
WRITE(6,1000 )     ISTAGE, (KPLANE( I )rIPLANE( I ), T=1,M) 

DP    STAGES    2    THRU   N 

DO    900    ISTAGE=2,N 

READ    IN    ARC    PARAMETERS    (REMAINING    ARCS) 

READ(5,60)    B,W,BLO 
BLOSUM=BLOSUM+BLG 
IF(BLOSUM.GE.UPBND)    GO   TO    2550 

SET    STATE    VARIABLE    &    RUN    FIB   SEARCH 

DO  800  J  =  2.M 
IX=J-1 

DO  200  N0=l,20 

IFdFIBNC(NO).GE.J)  GO  TO  300 
200  CONTINUE 
300  NCM=N0-1 

IA=0 

IB=IFIBN0(N0)-1 

DO   500    ITER=1, NOM 

N01=N0-ITER 

N02=N0-ITER-1 

IF(IA  +  IFIBNC(N01).GE.J  )    GO    TO    400 

ARG1=-B* (IA+IFIBN0(N01 ) J 

IT1=IX-IA-IFIBN0(N01)+1 

IF(N02.EQ.O)    GO    TO    350 

ARG2=-B* (I A+IFIBN0(N02 ) ) 

IT2=IX-IA-IFIBN0(N02)+1 

GO    TO    360 
350   ARG2=-R*IA 

IT2=IX-IA+1 
360    F1=W*£XP(ARG1)+FCN(IT1) 

F2=W*EXP(ARG2)+FCN(IT2) 

IF(F2.LT.F1)    GO    TO    400 

IA=IFIBN0(NO2)+l+IA 

GO  TO  500 
400  IB=IFIBN0(N01)-1+IA 
500  CONTINUE 

FCNNEW(J)  =  AMINKF2,F1) 

IPLAME( J)=IA 

IFUSTAGE.ZQ.N)  CAP (J) =FCNNEW ( J) +BLOSUM 

IFUSTAGE.EQ.N.AND.J.GT.2) 
1DELCAP( J)=FCNNEW(J)-FCNNEW(IX) 
800  CONTINUE 

FCNd  )=FCN(1  )  +  W 

DO    850    1=2, M 

FCN(I)=FCNNEW( I) 
850    CONTINUE 

WRITE (6,1000)    ISTAGE,(KPLANE(I),IPLANE(I) ,1=1, Ml 
900    CONTINUE 

CAPd  )=FCN(1  )+BLOSUM 

DELCAP(2)  =  FCN(2)-FCN(1) 


END  DP  STAGES 
1000  FORMAT( • !•////////"  f.29X.» STAGE  NUMBER  ».I2///// 


•  ,14X,3  CIST  ATE',  4X.  « IDECt ,  2X  )/// 


- ,14X,3( I6,4X,I4,2X) /• 
•,14X, 3( I6.4X.I4.2X)/' 

•  ,14X,3( 16, 4X, 14, 2X )/' 

•  ,14X,3( I 6,4X,I4,2X)  /» 
•,14X,3< 16, 4X, I4,2X)/» 

•  ,14X,3(I6,4X,l4,2X)/» 
•.  14X.3U6.4X.I4.2X)  /• 


-  ,14X,3(I6,4X,I4,2X)/ 
■  »14X,3(I6,4X,I4,2X)/ 
•tl4X, 3( 16, 4X, 14, 2X)/ 
• ,14X,3(I6.4X.I4.2X)/ 
•tl4X,  3(  I6,4X,I4,2X)/ 
1  ,14X,3(I6,4X,I4,2X  )/ 
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